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Why does the apparent mass of a coronal mass ejection increase? 

Li Feng^’^, Yuming Wang^’*, Fang Shen^, Chenglong Shen^, Bernd Inhester^, Lei Lu^, 

Weiqun Gan^ 

ABSTRACT 

Mass is one of the most fundamental parameters characterizing the dynamics 
of a coronal mass ejection (CME). It has been found that CME apparent mass 
measured from the brightness enhancement in coronagraph images shows an in¬ 
creasing trend during its evolution in the corona. However, the physics behind it 
is not clear. Does the apparent mass gain come from the mass outflow from the 
dimming regions in the low corona, or from the pileup of the solar wind plasma 
around the CME when it propagates outwards from the Sun? We analyzed the 
mass evolution of six CME events. Their mass can increase by a factor of 1.6 
to 3.2 from 4 to 15 Rs in the held of view (FOV) of the coronagraph on board 
the Solar Terrestrial Relations Observatory (STEREO). Over the distance about 
7 to 15 Rs, where the coronagraph occulting effect can be negligible, the mass 
can increase by a factor of 1.3 to 1.7. We adopted the ‘snow-plough’ model to 
calculate the mass contribution of the piled-up solar wind in the height range 
from about 7 to 15 Rs. For 2/3 of the events, the solar wind pileup is not suffi¬ 
cient to explain the measured mass increase. In the height range from about 7 to 
15 Rs, the ratio of the modelled to the measured mass increase is roughly larger 
than 0.55. Although the ratios are believed to be overestimated, the result gives 
evidence that the solar wind pileup probably makes a non-negligible contribution 
to the mass increase. It is not clear yet whether the solar wind pileup is a major 
contributor to the hnal mass derived from coronagraph observations. However, 
our study suggests that the solar wind pileup plays increasingly important role 
in the mass increase as a CME moves further away from the Sun. 
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Subject headings: Sunicoronal mass ejections (CMEs) — solar wind — Sun: 
evolution 


Introduction 


Mass is a major parameter characterizing the dynamics of a CME. However, due to the 
limitation of observations, we do not know very much about the CME mass, and its evolution 
from the corona to interplanet ary space. CME mass es are mostly derived from images of 


white-light coronagraphs (e.g., IVourlidas et all I 2 OOOI) . as a CME can cause an enhanced 


brightness in white-light images. However, the mass derived from white-light coronagraph 
images is just an ‘apparent’ mass of a CME. It has two components: one is the mass of 
plasma really ejected from the corona, the other is the piled-up solar wind mass due to the 
compression of the solar wind plasma surrounding the CME. It is difficult to separate the 
two components from imaging data. Therefore, what we often calculate is the apparent 
mass. A few studies have shown that the appa rent mass of a CME increases with time (e.g., 
Colaninno fc Vourlidas 2009 : Bein et ah 2013 ). 


The statistical work by IWang et ahl (1201 ll) has shown that the CME brightness is pos¬ 
itively correlated with the CME speed. A natural explanation is that a faster CME causes 
higher compression of ambient solar wind and therefore a brighter signature in coronagraph 
images. Numerical simulations also revealed that the mass of the high-density shell ahead 
of an erupted flux rope increases as it propagates away from the Sun b ut the mass of the 


dark cavity, whe re the flux rope is located, remains almost unchanged flLngaz et ahl 12005 


Rilev et al.N2008[) . All these studies imply that the piled-up solar wind mass may not be 
negligible, and may contribute signihcantly to the aforementioned phenomenon — apparent 
mass gain of a CME during its propagation. 

On the other hand, there is also evidence suggesting that t he CME mass gain can be 


attributed to the continuous outfl o w from the low c orona (e.g., iHoward fc VourlidasI 12005 


Bein et al.ll2013l: iTian et al.ll2012 ). iBein et ahl ( 20131) investigated the CME mass evolution 
for a set of 25 events, and found that the apparent mass incre ased by about 2% t o 6% per solar 


radiu s based on the COR2 observations taken by STEREO flKaiser et al.ll2008l: iHoward et ah 


20081 ). They further found that the mass centroids of these CMEs moved toward the Sun 
relative to their self-similar expansion in most events. The authors suggested that the reason 
might be that the rear part of the CME with higher mass might propagate slower than the 
front part of the CME. They found that the higher mass at the rear was s upplied by the 
outflow from the lower corona. However, in the analysis of iBein et ahl (120131) . the fact that 
materials may fall back to the Sun, which can also cause a backward motion of the CME 
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centroid, was not considered. Thus, the question how much pileup and outflow contribute 
to the mass increase still remains unsolved. 


In this paper, we are going t o answer the quest ion why the CME apparent mass increases 
in a way different from that bv iBein et al.l (120131) . We try to quantitatively figure out at 
which height and how significantly the piled-up solar wind mass contributes to the CME 
apparent mass. In Section O we briefly describe the selection of events, their observations, 
and their 3D geometric and kinematic properties. How we derive the CME mass and the 
mass evolution of selected events is presented in Section [31 In Section jH we analyze the 
possible source of the CME mass with the help of a solar wind pileup model and MHD 
numerical simulations. The final section consists of the discussions and conclusions. 


Events 


2.1. Event Selection 


Six CM E events ('se e Table fH) were selected from the halo CME list (http: //space . ustc. edu. cn/drea 
compiled bv IShen et al.l (l2013al) . The criteria to choose these events are: (1) They appeared 
halo in Large Angle Spectrometric Coronagraph (LASCO) on board Solar and Heliospheric 
Observatory (SOHO) spacecraft. They were also observed by CORl and COR2 corona- 
graphs on board STEREO and in most cases as almost limb events. Limb CMEs are less 
occulted by the coronagraph occulter, thus are favorable for mass estimate. (2) There was 
no contamination from other CMEs in the region where the investigated CME was located. 

(3) There was no contamination of very bright Ha line emission from prominence material 
inside the investigated CME in the coronagraph images. Mass estimat e of such contami¬ 


nated CMEs is oft en misleading if Thomson scattering theory is applied (IMierla et al.ll2011 


Carlev et aLll2012h . The six CMEs cover a broad range of leading-edge speed from 450 to 
1300 km s“^ as indicated in Table [U 


The twin spacecraft STEREO A-l-B and SOHO monitor the solar corona from three 
different angles of views at about 1 AU. CORl and CO R2 on board STERE O observe the 
corona from 1.4 to 4 Rs and 2.5 to 15 Rs, respectively ( Howard et ah 2008 ). LASCO C2 
and C3 cameras on bo ard SOHO have a FOV of 2 to 6 Rs, and 3.7 to 32 Rs, respectively 
(iBrueckner et al.lll995l) . 
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2.2. Deprojected CME Geometry and Kinematics 


The projection of a CME onto an image may significantly distort its geometric and kine¬ 
matic parameters. For a better estimate of the mass of a CME and its evolution, we have 
tried to derive the deprojected three-dimensional (3D) geometry and kinematics. Various 
models/methods have been developed to retri eve realistic parameters based on multi-p oint 


imaging data, e.g., the forward modeling fe.g.. Thernisien 

et al. 

2009 

Wood et al. 

2009). tri- 

angulation metho 

d fe.g., Liewer et al. 200‘ 

h Temmer et al 

. 200£ 

!), po 

larization ratio method 
letric localization (e.g.. 

(e.g., 

Moran et al 

. 201C 

i), inverse method ( 

Antunes et al. 

ZOOol). geon 

de Koning et ah ^ 

>009; 

3vrne et al. 2010: F 

eng et a 

2012) 

. The 

comparison of various mod- 

els could be found in. 

e.g.. Mierla et al.l 

(2010). 

riiernisien et al. 

(2011) and 

Feng et al. 


(l2013a|). 


Here we use the forward modeling developed bv lThernisien et ahl (120091) to derive the 3D 
geometric and kinematic parameters of the six CME events. In this method, the graduated 
cylindrical shell (CCS) model is applied, which consists of a tubular section forming the 
main body of the structure attached to two cones that correspond to the ‘legs’ of the CME. 
This model has six free parameters: the propagation longitude and latitude (0 and 9), aspect 
ratio {k), tilt angle( 7 ), half angle between two legs (a), and heliocentric height of the CME 
leading edge {h). From the parameters a and k, we further derived the edge-on angular 
width (2(5) and face-on angular width (2(q; -|- (5)), where 6 = arcsin(fi;) is the half angle of 
the cone, i.e., the CME leg. Note that all the heights mentioned in this work are the radial 
distance measured from the solar center. These six parameters were tuned to obtain a best 
match between the modeled flux rope and the observed CME from different viewpoints. More 


detail s of this model and its ap plica tion to CMEs co uld be found in, e.g., iThernisien et ah 


(120091) . IThernisien et al.l (120111) and iGui et ahl (120111) . For the time sequence of images of 
each CME, we obtained a set of parameters at each time. For simplicity, the free parameters 
are set to be constants by trial and error, except the height h which is a function of time. 


As a demonstration. Figured] shows the results for the event on 3 April 2010. The upper 
panels are the images taken by COR B, C2 and COR A, respectively, in which a pre-event 
image has been subtracted to remove the background and make the CME structure more 
pronounced. To get a picture of the CME as complete as possible, we incorporated CORl 
and COR2 observations. The aforementioned procedure has been applied to all the other 
Eve CME events as well. The lower panels are the same as the upper panels except that 
the green meshes are superimposed to show the best-ht flux rope. It has to be pointed out 
that for the CCS reconstruction, we did not include in the htting process the shock-sheath 
region when visible. By repeatedly applying the CCS model to the time sequence of the 
images of the CME, the evolution of the height of the leading edge of this CME can be 
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obtained and is shown in Figure |2j The uncertainty in the height is indicated by the error 
bars. It has two sources; one is the ^ 0.3 Rs uncertai nty in the definition of the CME leading 


front in the white-light images fe.g.. lBein et al.ll2013f) . and the other comes from the forward 
fitting process. We assumed that the uncertainty is about 5% of the height at each time 
instance. At the height of 20 Rs, it produces an uncertainty of 1 Rs wh ich is close to the 
maximal uncertainty in the sensitivity analysis in iThernisien et al.l (120091) . The height-time 
data points were further fitted by a linear function to estimate the average speed of the 
CME. 


Table [T] lists the geometric and kinematic parameters of our six events, i.e., longitude 
and latitude of the propagation direction (0 and 9), face-on and edge-on angular width (IT/ 
and We), CME solid angle (hi), and average speed (vavg). The longitude and the latitude are 
in a heliographic coordinate system with Earth sitting at zero longitude. Positive/negative 
longitude represents a CME propagating westwards/eastwards. The solid angle spanned by 
a CME is approximated by the multiplication of the face-on and edge-on angular width. The 
leading-edge speeds of the CME events vary from about 450 to 1300 km s“^, showing good 
representativeness. 


3. Evolution of CME Apparent Mass 
3.1. Mass Estimate Method 


Mass calculations base d on white-light coronagraph images are usually based on the 


(Minnaertlll930: 

van de Hulstlll950: 

Billines 

1966; 

Howard & Tannin 


20091 ). The major non-Thomson scattered contributions like stray light and dust scatter are 


removed by subtracting a pre-event background image. The remaining image brightness then 
should show the CME density distribution produced by Thomson scattering alone. For the 
total brightness in a given image pixel, the details of this scattering process can be described 
by 


b = bQ^ f F{r,x,u) Ne{r) de, (1) 

^ Jlos 

(1 - u)C{r) + uD{r) {1 - u)A{r) + uB{r) . 2 

F(r, Y, u) = 2- -^^ 1^ sm^ v 

^ ^ 1 - m/3 1 - m/3 ^ 

where 6© is the physical mean solar brightness (MSB) that is used as a unit of the pixel 
value in calibrated CORl and COR2 observations, F{r,x,u) is a function of three variables 
r, X, and u. r is the distance of the scattering location to the solar center, y is the scattering 
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angle between the radial vector from the solar center to the scattering electron and a position 
vector from the electron toward the observer, u accounts for the solar disk limb darkening, 
and we used the value of 0.56 for u in the visible white-light spectral range. A, B, C and D 
are known functions of r. is the differential Thomson scattering cross section, and has a 
value of 7.95 x 10“^® cm^ sr“^. 


In Equation ([T]), particular attention should be paid to the scattering angle x- In this 
work, the scattering angles were computed under the assumption that a CME is concentrated 


as in IColaninno fc Vourlidasl fl2009l) and ICarlev et al.l (120121) . And the CME propagation 
direction was derived in Section 2.2 using CCS forward modelling method. Under such all-in- 
propagation-plane assumption of a CME, we can calculate the column densi ty Np dl for 


any L OS with Equation ([T]). Assuming a composition of 90% H and 10% He fjVouriMas et ah 


20001 ). the mass in each pixel of a coronagraph image can be derived by multiplying the 


column density with the pixel area measured at a distance of about lAU from the observer. 


In our practical procedure, we hrst dehned a sector that contains the whole CME. In 
Figure El the mass in each pixel within the sectors in COR A and COR B is presented for 
the CME occurred on 3 April 2010. And we integrated the mass over the entire sector and 
obtained two mass values from COR A and COR B data, respectively. It should be noted 
that as we subtracted a pre-event background from the image of interest, the brightness in 
the non-CME region is supposed to be roughly at noise level. Thus, although the sector is 
larger than the region that the CME actually occupies, the calculated values of the apparent 
mass should be acceptable. Due to the different perspectives of COR A and B, the two 
values are slightly different. In Figure |4] we show the mass evolution with time. The mass 
evolution derived from COR A observations is delineated by the red solid line, and from COR 
B observations by the red dashed line. The average of the mass derived from COR A and B 
observations is delineated by the green solid line. We found that the difference between the 
mass derived from COR A and B was about 10 ±4%. Some other curves and associated error 
bars in Figure 0] will be explained in Section 3.2. The evolution of the averaged mass derived 
from COR A and B data for all six events is shown in Figure |5l The mass uncertainties will 
also be described in Section 3.2. 


3.2. Mass Uncertainty Analyses 


Vourlidas et al.l (120101) has shown that the instrumentally photometric uncertainty is 


quite small, and the subtraction of a pre-event ir nage introd u ces in signihcant errors, if the 
pre-event image is selected carefully. Similar to iFeng et al.l (120151) . we mainly considered 
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two sources of error in the mass estimate. One is the LOS distribution of the CME mass 
in Equation ([1]), especially the longitudinal propagation angle 0. The other source is our 
somewhat arbitrary selection of the CME sector. To estimate the uncertainties from these 
two sources, we varied 0 within 0 ± 5° at seven values and repeated the CME sector selection 
process independently for eight times. A total of 15 measurements were made from COR 
A and B observations, respectively. They are shown by black solid and black dashed lines 
in Figure 01 Correspondingly, the red solid and dashed lines are the averaged values over 
15 measurements. For our analyses of the mass evolution of a CME event, we adopted the 
green curve. 15 measurements from COR A and 15 measurements from COR B yield 15 
averaged mass over A & B. The green line and associated error bars are the mean and 3a 
values of the 15 averaged mass measurements. 


Figure 0] presents an example of the mass evolution and error analyses. In Figure O the 
mass evolution and associated error bars of all six events are shown. The black error bars are 
the 3a of all 15 mass measurements with different sectors and propagation longitudes, while 
the red error bars are the 3a of the mass measurements only with different sector regions. 
In most cases the major uncertainty comes from the selection of the CME sector. However, 
for the CME observed on 23 February 2012, the major uncertainty originates from the error 
of the propagation longitude. 


To explain the results of uncertainty above, we include in Figure [6] the dependence 
of the function F{r,x,u) on the scattering angle y at three impact distances of 2, 4, and 
10 Rs. The ordinate is F{r,x,u) normalized by its maximum at the PCS y = 90°. For 
COR observations, the impact distance is the projected distance of r from the scattering 
location to the solar center onto the POS. As revealed by Figure [6l when y deviates from 
90° within about 30°, F(r, y, u) only cha nges slightly and forins a pk teau around y = 90°. 
This phenomenon has also been found by Howard fc DeForest ( 201^ where they te r med i t 
“Thomson plateau” for the heliospheric HI observations and IColaninno fc VourlidasI (120091) . 
Therefore, a small change in propagation longitude 0 ± 5°, which causes a small change of 
y ± 5°, does not effect F{r,x,u) much within the 30°-deviation range. Subsequently, the 
column density and mass for a given pixel derived from Equation ([T]) does not change much 
when y deviates from 90° within about 30°. This is the case for the hrst hve CME events. 
Therefore, the uncertainty of the mass for these events mainly comes from the selection of 
the CME sector rather than the propagation longitude. 


In the case of the CME occurred on 23 February 2012, the propagation direction deviated 
from the POS as seen by STEREO A by about 50°, and from the POS as seen by STEREO B 
by about 90°. As the CME appeared in the FOV of COR B as a full halo, the corresponding 
mass calculation subjects to a large error. Therefore, we only used COR A data for mass 















analyses. When we derived the mass from COR A data, a small change in y yielded a large 
change in mass. It is because the slope of F{r,x,u) at y ~ 140° is quite steep in Figure |6l 
a small change of y ± 5° makes a large change in F{r,x,u), hence in column density and 
mass. That is why the major uncertainty of this event comes from the uncertainty of the 
propagation longitude. 


Note that our all-in-propagation-plane assumption yields an inhnitesimal angular width 


for a C ME. Mass can be underestimated under such an assumption according to lVourlidas et ah 


(I2OOOI) . The authors derived the ratio of the mass estimate under the infinitesimal-angular- 
width assumption to the mass under a finite-angular-width assumption as a function of the 
angular width of a CME. As the angular width usually keeps nearly constant in time in the 
COR FOV, the under-estimate ratio should also keep constant. Therefore, the mass change 
rate dm/dt is not effected by the inhnitesimal-angular-width assumption. 


3.3. Treatment of the Occulting Effect 


Figure [5] clearly shows that the mass evolution has an increasing trend. iBein et ah 


(120131) have demonstrated that one of the major causes of such an increase is the entering 
of the CME ejecta into the COR FOV from behind the occulter, which we call occulting 
effect. They then showed that the apparent mass evolution could be fitted by the following 
function 

hoi 


m{h) = mo 


1 - 


-h Am(h - hocc), 


( 2 ) 


where the parameter hocc is the height of the outer edge of an effective occulter, mo is a 
modeled initial apparent mass of a CME, and Am is an additional mass increment per unit 
height. The first term on the right-hand side of Equation ([2]) represents the increase of visible 
mass from a homogeneous cone-shaped and self-similarly expanding CME. It describes how 
the visible mass of a CME with a constant mass mo evolves, when entering the FOV of COR 
due to its expansion. The increase of the visible mass above the occulter is merely a geometric 
effect. The second term is an empirical correction to this simplihed model evolution, as the 
CME mass could increase due to phy sical reason s , e.g., mass outflow from a dimming region. 


or solar wind pileup. In the model of iBein et al.l (120131) . the actual physical mass increase of 
a CME can be approximated by m{h) = mo -|- Am(h — hocc)- 


In Figure [TJ the black discrete data points are the calculated mass as a function of 
deprojected height with the error bars in horizontal and vertical directions. Taking these 
error bars into account, we fitted Equation ([2]) to the mass-height measurements of each CME 
event. The best-ht parameters for each CME event are written in the bottom right, and for 
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mo and Am, the mean ±a are in logarithmic scale. In each panel the htting results are shown 
by the black solid line, and the black dashed line is the mass variation mo + Am{h — hocc) with 
the occulting effect removed. The vertical line marks the height from which the occulting 
effect starts to be negligible. It is dehned as the height from which the measured mass can 
reach above 97% of the dashed line. It is found that the occulting effect could be ignored 
once the CME reached about 7 to 8 Rs. We denote this height as hi, the corresponding 
mass at hi is mi. mi and hi for all six CME events are listed in Table |2l The other two 
mass quantities that are also included in Table |2] are mg below hgcc, and the hnal apparent 
mass m2 before a CME leaves COR FOV at h 2 . 

Note that our purpose of adopting this model is only to obtain the height where the 
occulting effect starts to be negligible. We did not use this simplihed model for further 
analyses of mass evolution. The measured mass evolution may deviate from the linear 
increase assumption in Equation (j2]). Actually, in Figure [71 for most events, when the CME 
starts to leave the COR FOV, the measured mass is lower than this linear-model prediction. 


4. Sources of CME Apparent Mass Gain 


As mentioned before, the CME apparent mass might consist a signihcant contribution 
from ambient compressed solar wind plasma. In some appropriate circumstances, shock 
might form around a CME moving faster than the fast magnetoacoustic wave. For events 
E3, E5, in COR images we saw deflected streamers which were not adjacent to the CME, 
and left a very diffusive s pace between the deflected streamer and the CME . It might be 
the signature of a shock (jVourlidas et al.ll2003l: [Ontiveros fc Vourlidasll200911 . It is almost 
impossible to separate the brightness contribution of the shock from that of the CME. 
Although solar wind pileup and compression due to a MHD shock are different physical 
processes, as the observed shocks were far fainter than the CME, we did not remove these 
possible shocks from our mass analyses and treated them as an extreme case of compression. 


To quantify the piled-up mass of the solar wind plasma surrounding the CME, a highly 
ideal model, called the ‘snow plough’ model, is applied. In this model, a CME accretes mass 
as it propagates and its momentum is changed through the interacti on with the a mbient 
solar wind. The mass change per unit time is described by (see also in [TappinI [20061) 


dm 

dt 


Psw-^\'^ sw ^ 


(3) 


where m, v and A are the mass, speed and the cross-sectional area of a CME, respectively, 
and Vsw and psw are the ambient solar wind speed and density, respectively. The value of 
A is approximated by the solid angle hi in Table [T] times h?, i.e.,A = 4(5(a -|- (5)h^. The 
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derivation of the above equation could be found in Appendix It should be noted that the 
equation assumes that the entire piled-up plasma moves together with the CME. Although 
this assumption may not be true, it is still worth investigating if this model can represent 
the mass gain during the CME propagation. Actually the ‘snow plough’ equation gives the 
maximal amount of mass that the solar wind can contribute to CME mass. 


4.1. Parameters of Background Solar Wind 


The ‘snow plough’ model requires the solar wind parameters. Thus we use 3D MHD 
simulations to obtain the undisturbed solar wind parameters along the CME path from 2.5 
to 20 Rs. The solar wind density and speed at different heights are calculated as the average 
over a spherical sector of 60 degrees which is centered at the 3D propagation direction of 
a CME. The numerical scheme in the 3D MHD sin iulations is a total variation diminishing 
(TVp) scheme in a Sun-centered spherical system (IFeng et al.ll2003L 120051: IShen et al.l 120071. 


20091) . For all six events, the synoptic maps of longitudinal magnetic held over a Carrington 


rotation from Wilcox Solar Observatory are used as input to the code for reality. The details 
of how to 


Shen et ah 

2613b: 

Wang et ah 

2014) 


The averaged solar wind density and speed against height over the spherical sectors are 
presented in Figure [8] and [9l respectively. For all six events, the variation of density follows 
about h~^ to h~^ below 7 to 8 Rs, and follows above 7 to 8 Rs. It is also the height where 
the occulter effect starts to be negligible. The solar wind speeds have a rapid increase at 
lower heights and then a slower increase at larger heights. The speeds have a similar trend 
of evolution for different Carrington rotations, but reach different values at 20 Rs which 
range from about 220 to 460 km s“^. The solar wind parameters at 20 Rs have been listed 
in Table [2] for reference. We also plot in each panel the speed of the CME derived from a 
quadratic £t to the height-time diagram. The scale of the solar wind speed is marked on the 
left Y axis, and the scale of the CME speed is marked on the right Y axis. For all events, 
we End that Ucme — "^sw decreases with height. 


4.2. How Much Mass Can Solar Wind Pileup Contribute? 

The calculation of the mass that the solar wind pileup can contribute is straightforward, 
as all the parameters m, n. A, and Usw can be obtained from either CME observations 
or MHD numerical simulations. The CME speed v was estimated from a quadratic £t to 
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the height-time plot, as it provides a better result than the linear £t. We started the mass 
pileup from rrii at hi, below which we assumed that the mass increase is due to the geometric 
occulting effect and the mass supply from lower corona. If the mass increase above hi was 
solely because of the solar wind pileup, we would obtain the mass evolution indicated by the 
red line in each panel of Figure [3 The shadow region in orange is obtained by assuming 
that the solar wind density and velocity have an uncertainty of 20%. 

It is found that, for most events (E3 through E6), the solar wind pileup mass is not 
sufficient to interpret the measured mass increase from COR observations represented by 
black symbols, suggesting that apparent mass increase is also contributed by the outffow 
from the low corona. For the event El, the solar wind pileup mass has a higher value than 
the measured ones above hi in the late-phase evolution in the COR FOV, while for the event 
E2, it is more or less consistent with the measured mass. Note that the ‘snow plough’ model 
gives an upper limit of the mass estimate. This model assumes that the solar wind plasma 
colliding with the CME will all be attached to the CME. However, this is an ideal case. 
Some solar wind plasma may not be attached, and just flow around the CME and eventually 
become the ambient background. Therefore, the actual mass increase due to the solar wind 
pileup may drop below the red lines. 

Although the solar wind pileup mass may be overestimated, the results still implies 
that the solar wind pileup probably makes a signihcant contribution to the apparent mass 
increase observed by the coronagraphs. By comparing the piled-up mass rris from hi to h 2 
with the total mass gain m 2 — mi during the same period, we End that the piled-up mass 
occupies more than half of the total mass gain as listed in the column for fs in Table |2l or 
ms is at least 19% of the hnal mass m 2 before the CMEs leave the COR FOV. Since the 
piled-up mass are overestimated and our inhnitesimal-angular-width assumption may make 
the measured mass underestimated, the ratios of the piled-up mass to the total mass gain 
given above are obviously overestimated. However, even if we consider a 100% overestimate, 
the ratios of the piled-up mass is larger than 25% in the height range from hi to h 2 . That 
is to say, although the pileup is not sufficient to explain the observed mass increase, its 
contribution is non-negligible. Due to the overestimation, it is not clear whether the pileup 
could make a major contribution to the measured mass increase. 

Besides, from Figure [71 we can End that the slopes of mass evolution dm/dh of the 
measured mass and modeled pileup mass gradually decrease with height. When the CME 
was about to leave the COR FOV, the two slopes are nearly the same for most of the events 
except the Erst one. It indicates that solar wind pileup makes a more important contribution 
to the mass increase at larger distances from the Sun. 
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5. Conclusions and Discussions 


Based on the brightness enhancement in coronagraph images, we have followed the mass 
evolution of six CMEs. The deprojected speed of the investigated CMEs covers from 450 
to 1300 km s“^. All these CMEs with different kinematics have an increasing trend in their 
mass evolution even when the occulting effect is removed. The long-lasting accumulation of 
CME mass in coronagraph FOV implies that the initial kinetic energy in the CME source 


region is smaller than the energy estimate using the hnal asymptotic mass (jCarlev et ah 
2012 : Feng et ah 2013b). 


Physically, there may be two sources of CME mass gain in the corona: the solar wind 


(Jin et ah 

2009; 

Tian et al. 

2012: 


We 

calculated the mass contribution of the solar wind pileup from the height beyond which 
the occulting effect is negligible. It is found that solar wind pileup probably may make a 
non-negligible contribution to the apparent mass increase observed by coronagraphs. For all 
of the events, the solar wind piled-up mass might occupy more than half of the total mass 
gain during the same period from 7 or 8 Rs to the edge of the FOV of COR2. It has to 
be pointed out that those ratios represent upper limits of the pileup contribution. However, 
even if we consider a 100% overestimation, the ratio of the pileup mass to the total mass 
gain could be larger than one fourth in the height range from 7 to 15 Rs. We also hnd 
that the contribution from the solar wind pileup becomes increasingly signihcant as a CME 
propagates to larger distances from the Sun. In short. Our work reveals that the solar wind 
pileup is not sufficient to explain the measured mass increase in the COR FOV. However, 
its contribution is non-negligible. Due to the overestimation, it is not clear yet whether the 
pileup could make a major contribution to the measured mass increase. 


Whether the solar wind pileup is a major contributor to the hnal apparent mass in the 
COR FOV is not clear either. If we assume that the mass below 7 or 8 Rs came from the 
mass outhow from the dimming region, the ratio of ms/m2 reveals that the solar wind pileup 
comprises less than half of the hnal mass for all events. In this case, it may not be a dominant 
contributor to the CME apparent mass in the COR FOV. Whether this assumption is valid 
or not requires further investigation. Future work involving the mass estimate from the 
dimming regions using SDO/AIA multi-wavelength data will be pursued. If the solar wind 
pileup could also contribute to the mass below about 7 Rs, for some of the CMEs, the pileup 
may be a dominant contributor. 


Both the piled-up mas s derived from the ‘snow-plough’ model and the ‘virtual mass’ 
introduced in ICargilll (120041) are related to the solar wind mass density. However, they are 
theoretically and observationally diherent. From hydrodynamic point of view, the ‘virtual 
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mass’ is the ‘added mass’ of an accelerating body moving in a fluid, and is req uired to obtain 
the correct accelerating force in the momentum equation. The ‘virtual mass’ in ICargilll (120041) 
was derived under the assumption of an incompressible fluid, which means that the density 
in the surrounding solar wind does not change over time. The total brightness observed 
with a coronagraph is proportional to the density, as there is no density change caused 
by the added mass, we will not be able to detect it from base-difference images. On the 
other hand, the piled-up mass in the ‘snow-plough’ model is a compressional effect, and the 
resulting density enhancement could in principle be detected in the base-difference images. 
Furthermore, pileup occurs even when a CME is moving steadily without any acceleration. 


DeForest et al.l (120131) tracked the CME mass from STEREO/COY12 FOV to HI2 FOV, 


and found that the mass of the CME under investigation continuously increased until to 
the distance of about 0.7 AU, but the mass increase dropped with distance. It is consistent 
with the trend of evolution predicted from our mass evolution in the COR FOV. When 
the CME started to leave the COR FOV, the measured mass increase per height {dm/dh) 
seems to agree with the prediction of the ‘snow plough’ model. As an outlook of the work in 
this paper, we will test whether this model can explain the mass evolution in the HI FOV. 
Furthermore, if we assume a proper shape of CME at 1 AU, the mass residing in the CME 
and in the shock sheath can be approximately calculated from in-situ measurements, which 
can add another data point in the mass evolution prohle. 


SOHO and STEREO are projects of international cooperation between ESA and NASA. 
The SECCHI data are produced by an international consortium of NRL, LMSAL, and NASA 
GSFC (USA), RAL, and U. Birmingham (UK), MPS (Germany), GSL (Belgium), IOTA, and 
IAS (France). SDO is a mission of NASA’s Living With a Star Program. L.F. and W.Q.G. 
are supported by the grants from MOSTG (2011GB811402), NSFG (11473070, 11427803, 
11233008 and 11273065), NSF of Jiangsu Province (BK2012889), GAS(XDA04076101). Y.W. 
and G.S. are supported by the grants from MOSTG (2011GB811403), GAS (KZZD-EW-01- 
4) and NSFG (41131065, 41121003 and 41274173). F.S. is supported by the grants from 
MOSTG (2012GB825601) and NSFG (41174150 and 41474152). L.F. also acknowledges the 
Youth Innovation Promotion Association, GAS, for hnancial support. The work of B.I. is 
supported by DLR contract 50 OG 1301. 


A. ‘Snow^ plough’ model for CME propagation 

‘Snow plough’ model assumes that all the solar wind plasma colliding with the GME 
will be attached to the GME. Based on this assumption and considering the ID problem. 
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the mass exchange between the CME and the solar wind is given by 

Am = -Arusw = PswA\vsw - v\At (Al) 


or 


dm 

dt 


dlTlsiu 

dt 


Psw-^l'^^sw V 


(A2) 


where m, v and A are the mass, speed and the area of the cross-section of the CME, respec¬ 
tively, and msw, Psw and Vgw are the mass, density and speed of the solar wind, respectively. 
The cross-section area A = 4:6{a + 6)h‘^, where 26 is the flux rope edge-on angular width, 
2{a + 6) is the face-on angular width, and h is the height of the leading edge. Using absolute 
value of {vsw — v) means the solar wind pileup could occur at either the frontside or backside 
of the CME. When the CME speed is higher than the solar wind speed, the solar wind piles 
up mainly at the frontside of the CME; When the CME speed is lower than the solar wind 
speed, the pileup mainly occurs at the backside. For the CME events in this paper, we have 
CME speed always higher than the solar wind speed. Therefore the pileup mainly occurred 
at the frontside. 


The mass exchange causes the loss of the momentum of the background solar wind, 
which is 

Apsw ATTlg^Vsw PswA\Vs^ v\vsi^At ^A3) 

or 


^Psw dfriguj 

dt dt 


PswA\Vs 


VVg 


(A4) 


Meanwhile, the conservation of momentum requires ^ i.e., 

dmv dm.,,, 


dt 

dm 


dt 

-V + ma = ■ 


dm. 


dt dt 

where a is the acceleration of the CME. Using Eq lA21 the above equation leads to 

dm 

ma = -^V^sw - v) 


On the other hand, combine Eq iA4l and | 

dm 
dt 

After Substituting Eq lA7l for ma in the above equation, it is obtained that 

dm 


we have 
-n + ma = PswAvsvi\vsw - v\ 


dt 


PswA\Vg,,jj v\ 


(A5) 

(A6) 

(A7) 

(AS) 

(A9) 


which is the equation we used to calculate the mass contribution from solar wind pileup. 
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Table 1: The geometric and kinematic parameters of six CME events. 


No. 

Date 


on 

'N, n 

We (°) 

fl (sr) 

Uavg (km s“^) 

El 

2010-02-12 

-2.5 

7.2 

68 

47 

1.0 

752 ± 31 

E2 

2010-04-03 

4.6 

-25.1 

77 

50 

1.2 

824 ± 29 

E3 

2011-04-17 

-177.1 

8.9 

72 

52 

1.1 

1312 ± 51 

E4 

2011-06-21 

-20.8 

9.5 

70 

70 

1.5 

1006 ± 38 

E5 

2011-08-03 

15.4 

16.4 

68 

41 

0.8 

1294 ± 56 

E6 

2012-02-23 

67.2 

24.4 

139 

82 

3.5 

445 ± 12 


Note: (j) and 6 are the longitude and latitude of the 3D CME propagation direction, respectively. Posi¬ 
tive/negative values of b mean westward/eastward propagation with respect to the Earth. Wf = 2{a + S) 
and We = 2S are the face-on and edge-on angular width, respectively. D = W/We = 4:S{a -I-15) approximates 
the solid angle of a CME. Uavg is the average speed of the CME from the linear fit to the observed height-time 
profile. 


Table 2: Derived parameters of CME mass at different heights, solar wind parameters and 
its mass contribution 


No. 

mo 

hocc 

mi 

hi 

m2 

h2 


'^SW 

rrig 

fs 


logio(g) 

(Rs) 

logio(g) 

(Rs) 

logio(g) 

(Rs) 

(cm“^) 

(km s“^) 

logio(g) 


El 

15.17 

2.83 

15.27 

7.74 

15.37 

14.2 

869 

263 

15.09 

< 2.55 

E2 

15.38 

2.34 

15.48 

7.76 

15.60 

14.3 

657 

376 

15.01 

< 1.06 

E3 

15.35 

1.82 

15.63 

7.01 

15.86 

15.7 

824 

243 

15.32 

< 0.70 

E4 

15.45 

2.69 

15.62 

7.97 

15.78 

15.8 

696 

355 

15.20 

< 0.85 

E5 

15.42 

2.62 

15.68 

8.46 

15.85 

15.6 

878 

226 

15.15 

< 0.62 

E6 

15.65 

2.79 

15.83 

7.51 

16.01 

15.3 

610 

338 

15.28 

< 0.55 


Note: mo is the modeled initial apparent mass of a CME below the effective occulter height hgcc- mi is the 
apparent mass at the height hi where the occulting effect begins to be negligible, m 2 is the final apparent 
mass before a CME leaves COR FOV, and ^2 is the corresponding height. Ugw and Vgw are the ambient 
solar wind density and speed at h = 20Rs. mg is the solar wind mass piled up around the CME from hi 
to ^ 2 . fs= uicasures the contribution of the solar wind pileup to the apparent mass increase in the 

height range from hi to / 12 . 
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Fig. 1.— Upper panels: the CME observed by STEREO CORl & COR2 and by 
FORO/LASCO C2 on 3 April 2010. Note that CORl and COR2 images are incorporated 
to have a comprehensive overview of the CME morphology. Lower panels: The same as the 
upper panel except that the modeled flux rope indicated by green meshes are superposed. 
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Start Time (03—Apr—10 09:26:40) 


Fig. 2.— Deprojected height as a function of time for the CME event on 3 April 2010. The 
asterisks are the measurements and the black solid line is a linear ht. The average speed 
derived from the linear £t is marked in the upper part. The error bars indicate the 5% 
uncertainty of the deprojected heights. 
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Fig. 3.— Mass images within the CME sector regions. The plus sign and a white circle 
indicate the solar center and solar disk, respectively. The mass color bar is shown on the 
right side. 


mass (LoglO g) 
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Fig. 4.— Mass evolution with time for the CME on 3 April 2010. The thin black solid 
lines are the mass derived from COR A total brightness images with different selections of 
CME sectors and varing propagational longitudes within 0 ± 5°. A total of 15 independent 
mass measurements were made for the error estimate. The mass along the thin black dashed 
lines are derived in a similar way but from COR B observations. The mean mass evolution 
derived from COR A and COR B are marked by solid and dashed lines in red, respectively. 
15 measurements from COR A and 15 measurements from COR B yield 15 averaged mass 
over A & B. The green line and associated error bars are the mean and 3a values of the 15 
averaged mass measurements. 
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Fig. 5.— Mass evolution of all six events. The black error bars are the 3a of all the 15 mass 
measurements with different sector regions and different propagation angles. The red error 
bars are the 3a of eight mass measurements only with different sector regions. For most 
cases, major error comes from the selection of the CME sector. While for the CME on 23 
February 2012, the major error comes from the uncertainty of the propagation angle. 
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Fig. 6.— Dependence of the function F{r,x,u) on the scattering angle x at three impact 
distances (R=2, 4, 10 Rs). For COR observations, the impact distance is the projected 
distance of r from the scattering location to the solar center onto the POS. The brightness 
is normalized by its maximal value at the POS {F{r,x = 90°,m)). The vertical dotted line 
marks the position of y = 90°. 
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Fig. 7.— Mass profiles as a function of the deprojected height. The best fit result of 
Equation ([2]) is indicated by a black solid line in each panel. The mass evolution after 
removing the occulting effect is shown by the black dashed line. The best fit parameters 
are marked in the bottom right in each panel. For tuq and Am, the meania is written in 
logarithmic scale. The vertical dotted line marks the height hi where the occulting effect 
starts to be negligible. The mass evolution due to the solar wind pileup above hi is shown by 
a red dashed line in each panel. The uncertainties indicated by the shadow region in orange 
are derived under the assumption that and Vsw vary in the range from 80% to 120% of 
the simulated values. 
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Fig. 8.— The solar wind number density against height derived from MHD simulations. 
The number density is the average of the density over a cone of 60 degrees centered at the 
CME propagation direction. Three dashed lines indicate the variation of density following 
and h”®, respectively, from upper to bottom. The vertical dotted lines mark the 
height at 7.5 Rs. 
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Fig. 9.— Solid lines: the solar wind radial speed against height derived from MHD simu¬ 
lations, which is the average of the speed over a cone of 60 degrees centered at the CME 
propagation direction. Plus signs: the CME speed. 














